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ABSTRACT 

Context. The chemistry of proto-planetary disks is thought to be dominated by two major processes: photodissociation near the disk 
surface, and depletion on dust grains in the disk mid-plane, resulting in a layered structure with molecules located in a warm layer 
above the disk mid-plane. 

Aims. We attempt here to confront this warm molecular layer model prediction with the distribution of two key molecules for disso- 
ciation processes: CN and HCN 

Methods. Using the IRAM Plateau de Bure interferometer, we obtained high spatial and spectral resolution images of the CN J=2-l 
and HCN J=l-0 lines in the disks surrounding the two T-Tauri DM Tau and LkCa 15 and the Herbig Ae MWC 480. Disk properties are 
derived assuming power law distributions. The hyperfine structure of the observed transitions allows us to constrain the line opacities 
and excitation temperatures. We compare the observational results with predictions from existing chemical models, and use a simple 
PDR model (without freeze-out of molecules on grains and surface chemistry) to illustrate dependencies on UV field strength, grain 
size and gas-to-dust ratio. We also evaluate the impact of Ly a radiation. 

Results. The temperature ordering follows the trend found from CO lines, with DM Tau being the coldest object and MWC 480 the 
warmest. Although CN indicates somewhat higher excitation temperatures than HCN, the derived values in the T Tauri disks are very 
low (8-10 K). They agree with results obtained from C2H, and are in contradiction with thermal and chemical model predictions. 
These very low temperatures, as well as geometrical constraints, suggest that substantial amounts of CN and HCN remain in the gas 
phase close to the disk mid-plane, and that this mid-plane is quite cold. The observed CN/HCN ratio (^ 5 - 10) is in better agreement 
with the existence of large grains, and possibly also a substantial contribution of Ly a radiation. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - individual: LkCal5, MWC 480, DM Tau, - 
Radio-lines: stars 

1 ■ Introduction directly illuminated by the stellar UV and dominated by photo- 

dissociation reactions. The lowest layer, near the disk plane, 
^ Studymg the physical and chemical structure of proto-planetary ^^^^ ^^^^^ ^^^^ molecules are expected to be depleted on 
O . disks is a prerequisite for our understanding of planet formation. ^^^^^^ Between those, the third layer, where the gas is lukewarm 
^ . While much theoretical work has been made on the structure ^^e chemistry dominated by gas-phase reactions. Such a sim- 
and evolution of these disks (Berg in et al. || 2007i and references ^-^^^^^ ^^^^^^^ challenged by new observations which re- 
^ therein), key parameters such as gas density and temperature, or yeal that larg e amo unts of cold CO are still in the gas phase 
k> molecular content and distribution remain poorly constrained, (jpartois et al J l20Q3l: iPietu et aLlEoOTh . contrary to predictions. 
^ Since H2, which represents 80 % of the gas mass, cannot be y^^^^^^ ^-^^^^ proposed to refurbish CO in the c old zone 
^ „ observed with a good spatial resolution, our knowledge of the ( Aikawa & Nomura 2006; Aikawa 2007; ISemenov et a l. 2006, 
- - molecular component relies on low abundance tracers such as ^^so Willacy et al 2006) 
CO ( Guilloteau & Dutrev 1998). Interferometric maps of the CO 

emission, at mm and sub-mm wavelengths, resolve the gaseous ^his scheme, molecules are displaced from the upper lay- 

disks, revealing their radial and vertical structure. Studies of the ^o the cold mid-plane and have the time to cool down before 
rare CO isotopologues, ^^CO an dC^^O, confir m the existence of ^^ickmg onto grams; the results, however, are strongly dependent 
a vertical temperature gradient ("Dartois et a^^ iPietu et al.l «^ the grain surfacere^tion eflaciencies and on the disk struc- 
[2007h . allowing first face to face confrontation between model ture.| Hersant et al. | (|2009D have also shown that photo-desorption 
predictions and observations in the outer parts of the disks. ^crease the amount of cold CO m the gas. 

This led to a global picture of a flared disk consisting of 3 A deeper understanding of the disk structure relies on the 
strata, or layers. The upper layer, close to the disk surface, is study of molecules probing diflTerent physical and chemical con- 

ditions. Apart from CO, only few molecules show strong enough 

^ Based on observations carried out with the IRAM Plateau de Bure mm/sub-mm emission to be d etected: HCQ+, N2H^, HCN, CN 
Interferometer IRAM is supported by INSU/CNRS (France), MPG HNC , CCH, CS and H9CO (IDutrev et al.1 [ l997l: iKastner et afl 
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(Germany), and IGN (Spain). ' 1 19971: iThi et alJ 120041: Ivan Zadelhoff et all l200lh . Even fewer 
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can b e mapped ( Aikawa et al .1120031: lOi et al ] |2004l:lDutrev et alJ 
120071: ISchreyeretaLl 12008"). Among those, HCN and CN are 
particularly interesting because they trace different conditions 
and their abundances are closely related: HCN may be photo- 
dissociated (e.g., by Ly a radiation from the star) to yield CN 
(iBergin et al .1120031) . One thus expects HCN to trace the UV- 
shielded lukewarm and cold layers of the disk and CN the up- 
permost layer. 

In this article, we present high angular resolution (3'' to 
1.5'0 observations of HCN J= 1-0 and CN J=2-l emission in 3 
protoplanetary disks. The sources belong to the nearby Taurus- 
Aurigae star forming region, but have moved outside their parent 
cloud, so that they can be observed without interference. Table[T] 
summarizes their stellar properties. 

- DMTau is a classical T-Tauri star of mass 0.5 MQ_sur- 



rounded by a disk in Keplerian rotati o n (iGuilloteau & Dutreyl 
119981: lPietuetal.|[2QQ7b . iDartois et al.l ^2003) have shown the 
existence of a vertical temperature gradient from high angular 
resolution observations of several transitions of ^^CO ^^CO and 
C^^O. Eight different molecular species have been detected in 
this source (Dutrey et al. 1997). 

- LkCa 15 is a 1 Mp T-Tauri s t ar surrounded by a disk with 
a central cavity dPietu et al.ll2006h . iPietu et al.1 (l2007l) found no 
evidence of a vertical temperature gradient from their CO multi- 
isotopologue data. The cavity radius derived from dust emission 
is - 45 AU but some CO gas remains down to 20 AU or so 
dPietu et al.ll2007h . 

- Finally, MWC480 is a Herbig Ae star of 1.8 Mp sur- 
rounded by a Keplerian disk (iMannings et al.|[T997l: ISimon etaP 
I2OOOI: iPietuet al.1 l2007h . CO isotopologues indicate a vertical 
temperatu re gradient. The disk continuum emission was re- 
solved by iPietu et aP (l2006h . suggesti ng even lower tempera- 
tures for the dust than from ^^CO, but IGuilloteau et al.l (l201lh 
show that alternative solutions with radially varying dust prop- 
erties are possible. 

Section 2 presents the observations. The method of analysis 
and the results are shown in Section 3 and Section 4 is dedicated 
to the chemical modeling. We discuss the implications of our 
study in Section 5 and then we conclude. 



2. Observations 

All the observations reported here were made with the IRAM 
Plateau de Bure interferometer (PdBI). 

The CN J=2^1 observations of DM Tau were performed in 
Feb., Apr., Sep. and Oct. 1997 with baselines ranging from 15 
m to 180 m. The spatial resolution (HPBW) was 1.7 x 1.2'' (PA 
122°) for "natural weighting". The HCN J=1^0 data were ob- 
tained in Nov. 1999, Jan. 2001 and Jan. 2002 with baselines from 
15 to 330 m. The HPBW was 3.7 x 2.4'' (PA 39°) using natural 
weighting, and 3.0 x 1.8" at PA 45° for "robust" weighting. 

The LkCal5 and MWC480 data were observed in HCN 
J=1^0 and CN J=2^1 data in a track sharing mode, to obtain 
an homogeneous calibration for both sources. Baselines ranged 
from 15 to 330 meters, with an emphasis on short baselines: nat- 
ural weighting yielded an HPBW of 2.2 x 1.3" (PA 56°) and 
robust weighting an HPBW of 1.2 x 0.7" (PA 40°) for CN, and 
proportionally larger by a factor 2.5 for HCN.^ 

We used the GILDASQ software package ('Pety"2005') for the 
data calibration and deconvolution, and as framework for the im- 
plementation our minimization method. The resulting velocity- 



^ See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS software. 



channel maps in the HCN J=1^0 and CN J=2^1 lines are 
shown in Fig.[2]and[T]) and the spectra integrated over the disk are 
presented in Fig.|4l|3l In MWC480, the HCN emission is weak, 
but significant at the Icr level (see Fig|5]). Channel maps for HCN 
in MWC 480 line are only presented in Appendix, which also 
includes a more complete set of channel maps for all sources, 
in particular those of the other observed CN hyperfine compo- 
nents. We caution that although the integrated spectra in Fig.|4l|3] 
were obtained from the Clean images, they underestimate the 
flux. This is not due to the lack of short spacings, but to the dif- 
ficulty of deconvolving channel per channel low signal-to-noise 
data when the synthesized beam has relatively high sidelobes. 
Furthermore, we stress that this imaging bias has no influence 
on our derived parameters, since all our analysis is based on a 
comparison between model visibilities and actual measurements 
in the UV plane. Additional channel maps for CN, and also HCN 
in MWC 480 reveal 

Calibration errors could directly afl'ect the derived tempera- 
ture and molecular column densities. To minimize the calibration 
errors, all our flux density scale is based on a comparison with 
MWC 349, a sour ce with precisely known spectral index of 0.6. 
In addition, as in iGuiUoteau et al.1 (l201 ih . MWC 480 is strong 
enough in continuum to provide an additional comparison; the 
final relative flux scale accuracy is about 5%. Our flux density 
scale here is identical to that used by Guilloteau et al. (2011). 
However, recent re-calibration of the MWC 349 flux density at 
mm wavelengths indicate that this overall scale may be too low 
by about 10 % (Krips etal. in prep.). 

The velocity was smoothed to 0.41 km.s"^ in the cases of 
MWC 480 and LkCal5 to produce the velocity-channel maps 
on Fig. [Hand to 0.206 km.s"^ in the case of DMTau. We used 
however full spectral resolution in the model fitting analysis 
described below: for CN(2-1), 0.206 km.s"^ in LkCal5 and 
MWC 480, and 0.103 km.s"^ in DMTau. For HCN(l-O), spec- 
tral resolutions were 0.52 km.s"^ in LkCal5, 0.40 km.s in 
MWC 480 and 0.26 km.s"^ in DMTau. The velocity scale is re- 
ferred to the strongest hyperfine components, i.e. 226.874745 
GHz for CN, and 88.631848 GHz for HCN 

Prior to this paper, CN and HCN emission has been al- 
ready detected in several T-Tauri and Herbig Ae stars (DM Tau, 
GGTau; LkCa 15, HD 163296, MWC 480 and V4046 Sgr), but 
only with single-dish telescopes, hence at low spatial resolution 
(Dutrey et al. 1997; Thi et al. 2004; Kastner et al. 2008), except 
for the most nearby star TW Hya dOi et al.ll200 8). A preliminary 
study of CN(2-1) and HCN(3-2) at ^ 2. 5" resolution has been 
recently published by lOberg et al.l (l2010h . The observations rel- 
evant to our sources are summarized in Tabled 



3. Method of analysis and results 

Whereas the velocity-channel maps of Fig. [2] and [TJ were pro- 
duced from the dirty images via Clean-based deconvolution, the 
disk properties were derived directly from the observed visibili- 
ties through a model-fitting. For this, we used the disk model and 
hypothesis described in Pietu et al. ( 2007). We assume circular 
symmetry. The surface density of molecules follows in radius r 
a power law of exponent /?; perpendicularly to the plane of the 
disk, the density of molecules follows a Gaussian of half- width 
at 1/e intensity (scale height) h(r). We further assume the popu- 
lation of each molecular species can be represented by a rotation 
temperature T which depends only on r and follows a power law 
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Source 


Right ascension 


Declination 


Spec, type 


Effective temp.(K) 


Stellar lum.(Lo) 


CO paper 


MWC480 


04:58:46.264 


29:50:36.86 


A4 


8460 


11.5 


1,2 


LkCa 15 


04:39:17.790 


22:21:03.34 


K5 


4350 


0.74 


1,2 


DMTau 


04:33:48.733 


18:10:09.89 


Ml 


3720 


0.25 


1,3 



Notes. Columns 2&3: J2000 coordinates deduced from the fit of the PdBI 1.3 mm continuum map. Errors on the astromet r y are of order < 0.05''. 
Columns 4, 5, and 6, the spectral type, effective temperature and the stella r luminosity are tho se given in ISimon et al.l (l2000h . Column 6, CO 
interferometric papers are: l ^Pietu et al.. C2007h . 2 jSimon et al.. (2000.) . 3 = lDartois etaP ( l2003h . 



Table 2. Integrated flux of CN and HCN lines 





DMTau 


Sy (Jy km s 
LkCal5 


MWC480 


HCN J=l-0 


0.4^ 






CN J=l-0 


1.4^ 






CN J=2-l 


8.774.7^ 


6.8^ 


3.3^ 


HCN J=3-2 


2.9^ 


5.5^ 


2.3^ 


CN J=3-2 




12^ 


5.3^ 


HCN J=4-3 




4.6^ 


< 1.3^ 



Notes . Integrated line intensities from ^ IRAM 30m observation 
(lDutrevet al.ll997h and ^ JCMT observation (Thi et al. 2004). The con- 
version factor between Jansky to Kelvin (main beam temperature) for 
the JCMT (beam 13.5") at 345GHz is 18.2 Jy/K. ' SUA data from 
lOberg etal.l(l20TQh . 

of exponent q. The disk is thus described by 

T(r) = ro(r/7?o)-^ (1) 
^(r) = ^oir/RoT' (2) 
h(r) = Ho(r/Ro)-\ with /z = -1.25 (3) 

plus a constant turbulent width SV and an inner and an outer 
radius. If the molecule is at LTE, T(r) is equal to the kinetic 
temperature. 

3.1. Minimization technique 

The minimization tec hnique used to der ive the disk parameters 
has been described by lPietu et alJ (l2007l) . The continuum emis- 
sion was subtracted in the UV plane prior to the analysis, rather 
than being fitted sim ultaneously . 

As explained by Pietu et all (12007). temperature and surface 
densities can be derived by the minimization procedure from a 
single transition, even without hyperfine (hf) structure, under the 
assumption of power law distributions, provided the spatial res- 
olution is high enough to partially resolve the optically thick 
core. This is because the shape of the radial dependence of the 
brightness changes when the line becomes optically thick, from 
Th(r) = T(r) in the optically thick region, to a function of sur- 
face density and temperature in the optically thin domain (e.g. 
Th(r) oc ll(r)/T(r) at high temperatures). Lines with hf compo- 
nents yield additional information that may even allow to test the 
power law hypothesis. 

Handling of the hf components in the modeling requires ad- 
ditional care. Under our assumption of common excitation tem- 
perature, the apparent hf intensity ratios should depend only on 
the line opacities. However, when comparing model and obser- 
vations, limited spectral resolution and discrete velocity sam- 



pling may introduce a bias. First, the spectral backend delivers 
the flux integrated over a channel width, and not just sampled 
at the channel frequency. Second, the separation between the hf 
components is not an integer multiple of the velocity sampling. 
In the case of HCN, the intrinsic line widths are comparable to 
the channel resolution; then the difference between sampled and 
integrated values, combined with different sampling for each hf 
component, may change the predicted hf intensity ratios, which 
in turn leads to incorrect opacities and to erroneous tempera- 
tures and surface densities. We ensured that no such bias is intro- 
duced by comparing the results obtained directly from the disk 
model by using the same sampling in frequency as in the real 
data to those derived with a much higher sampling, followed by 
a smoothing that simulates the channel widths and center fre- 
quencies. We found almost no differences between the two sets 
of results, except for a very small (0.02 km.s"^) increase in the 
line width. The hf line intensity ratios were almost identical. In 
the case of CN, the intrinsic spectral resolution, 0.05 km.s"\ is 
small compared to the line widths. 

In the minimization procedure, all HCN hf components were 
fitted together assuming a common excitation temperature. For 
CN, the hf components were used in the analysis in two ways: 
all together and the intrinsically strong and weak components 
separately. The 3 sets of results were then compared. 

3.2. Analysis of CN and HCN lines 

In a first step of the analysis, all geometric parameters of the disk 
were left free. The derived para meter values were i n complete 
agreement with those derived by iPietu et aO (l200 7l) from CO, 
but had a lower accuracy as the lines are weaker. Accordingly, 
we fixed the inclination and position angle of the disk, as well 
as the stellar mass and gas scale height using Pietu et al] (l2QQ7h 
values (see Table [3]). For LkCa 15, we fix the internal radius Rint 
to 38 AU, a compromise between the CO and dust inner radii 
from Pietu et al. (2006, 2007). In a second step only So and p 
(the value and exponent of the surface density law), Rquu the 
outer radius, Tq and q (the value and exponent of the temperature 
law) and SV, the local line width of each molecule, were fitted. 

In addition, as many chemical models predict all molecules 
should be frozen onto dust grains close to the disk mid-plane, 
we tested models in which the molecules are present only 
above the disk plane. We tested two different parameterizations. 
Parametrization (1) is a simple geometric constraint, in which 
molecules only exist above z(r) = adh(r), with (the "de- 
pletion scale height") being a free parameter, i.e. being the 
molecular abundance, 

X^(r,z) = for z < aMr) (4) 

However, the chemical model of Aikawa & Nomural (l2006h sug- 
gests that, in a given disk, molecular abundances depend essen- 
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Fig. 1. Channel maps of the CN 2-1 emission for the strongest 
group of three blended hyperfine components. Top: LkCa 15; 
the spatial resolution is 1.7 x 1.0" at PA 44°, contour spacing 
is 60 mJy/beam, or 0.81 K and 2 Act. Middle: DMTau; the spa- 
tial resolution is 1.7 x l.T' at PA 120°, contour spacing is 53 
mJy/beam, or 0.61 K, 2.0o-. Bottom: MWC 480; the spatial res- 
olution is 2.2x 1.35'' at PA 60°, contour spacing is 55 mJy/beam, 
6r 0.45 K and la. 



Fig. 2. Channel maps of the HCN 1-0 emission. Top: LkCa 15; 
the spatial resolution is 4.1 x 2.9'' at PA 51°, contour spacing 
is 16 mJy/beam, or 0.21 K and 2. Oct. Bottom: DMTau; the spa- 
tial resolution is 3.7 x 2.4'' at PA 39°, contour spacing is 12 
mJy/beam, or 0.21 K, 2.0o-. 



tially on the (molecular hydrogen) column density between the 
disk surface and the location in the disk. Thus, at large radii, 
molecules get closer to the disk plane than at smaller radii. As 
our measurements are most sensitive to the 150 - 400 AU re- 
gion, assuming a constant ad as above may be inappropriate. We 
thus also considered a second option, where molecules are only 
present in the upper disk down a constant depth in (molecular 
hydrogen) column density, Ej, the "depletion column density". 
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Velocity (km.s ) 




Fig. 3. Integrated spectra for CN 2-1. The relative positions and 
intensities of the hyperfine components are indicated in red. Top: 
MWC480; middle: LkCa 15; bottom: DMTau 

Table 3. Modeling: fixed parameters 



Parameters 


DMTau 


LkCal5 


MWC480 


RA.(^) 


65 


150 


58 


in 


-35 


52 


37 


V/,,(km.s ^) 


6.00 


6.28 


5.02 


Vioo(km.s"^) 


2.15 


2.90 


4.05 


Rint(AU) 


1 


38 


1 



Notes. The parameters are taken from lPietu etaP (l2006lJ2007h . 
i.e. 

foo 

which in the isothermal approximation yields 

X„{r,z) = for ^ ^"^"(l^) ^ 



(5) 



(6) 




Velocity (km.s ) 





-5 
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Velocity (km.s ) 


10 


15 




Velocity (krr 



Fig. 4. Integrated spectra for HCN 1-0. The relative positions and 
intensities of the hyperfine components are indicated in red. Top: 
MWC 480, with the best fit model (see Tab|4|) superimposed to 
illustrate the weak detection; middle: LkCa 15; bottom: DMTau 



erfc being the complementary error function. This second option 
is likely to be more realistic, as it avoids placing molecules far 
off' the plane at large distances. Note that ad = corresponds to 
I^d being large (equal to the H2 column density). 

3.3. Results 

The disks of the two T Tauri stars, DM Tau and LkCa 15, display 
very similar properties (see Table |4l). Both show very low CN 
and HCN rotation temperatures, and their strongest hf line com- 
ponents have moderate optical depth around 100 to 300 AU. The 
fit to the weakest, optically thinner, CN hf components yields a 
nearly fiat surface density distribution and the fit to the stronger 
hf components a much steeper one. As a consequence of our 
single power law hypothesis, the temperature distributions fol- 
low the opposite trend, i.e. the temperature decreases faster for 
increasing r for the optically thin lines. Such a behavior may 
indicate that the single power law is not an appropriate repre- 
sentation of the CN distribution and that the density falls down 
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Table 4. Results of analysis 



Source 


Molecule 




P 




6V 


T, 


q 








(10^2 ) 




(AU) 


(kms-^ 


(K) 






MWC480 


HCN 1-0 


1.1 +0.4 


2.4 + 0.4 


[550] 


0.3 + 0.2 


[30] 


[0] 


0.0 + 0.7 




CN 2-1 (strong) 


11 + 1 


2.1 +0.1 


540 + 40 


0.25 + 0.04 


33 + 6 


[0] 


0.0 + 0.3 




CN 2-1 (weak) 


11 + 1.2 


2.1+0.15 


550 + 70 


0.26 + 0.07 


33 + 6 


[0] 


- 




CN 2-1 (all) 


10.4 + 0.9 


2.1+0.1 


545 + 35 


0.25 + 0.04 


30 + 4 


[0] 


- 


LkCal5 


HCN 1-0 


10.6 + 1.5 


1.1 +0.2 


600 + 40 


0.20 + 0.03 


7.0 + 0.6 


0.55 + 0.25 


0.0 + 0.8 




CN 2-1 (strong) 


38 + 4 


2.3 + 0.2 


570 + 20 


0.18 + 0.02 


10.8 + 0.7 


0.20 + 0.07 


0.15 + 0.17 




CN 2-1 (weak) 


115 + 18 


0.2 + 0.2 


690 + 20 


0.26 + 0.03 


6.6 + 0.3 


0.80 + 0.07 






CN 2-1 (all) 


58 + 5 


0.8 + 0.1 


630 + 35 


0.18 + 0.01 


8.8 + 0.3 


0.95 + 0.05 




DMTau 


HCN 1-0 


6.5 + 0.9 


1.0 + 0.3 


660 + 20 


0.18 + 0.01 


6.0 + 0.4 


0.00 + 0.12 


0.0 + 0.5 




CN 2-1 (strong) 


26 + 4 


2.1+0.15 


650 + 20 


0.16 + 0.01 


8.6 + 0.5 


0.05 + 0.05 


0.0 + 0.3 




CN 2-1 (weak) 


46 + 9 


0.7 + 0.15 


550 + 25 


0.20 + 0.01 


6.5 + 0.5 


0.50 + 0.08 


0.0 + 0.3 




CN 2-1 (all) 


35 + 9 


0.6 + 0.06 


620 + 15 


0.17 + 0.01 


7.5 + 0.3 


0.60 + 0.05 


0.0 + 0.3 



Notes. Column density 2(r) = 2 x {r/Ro)~P, external radius Rout, local line width 6 V, temperature T{r) - T x (r/ro)~^ and depletion scale height 
ad (see text). All values are referred to tq - 300 AU. Square brackets indicate fixed parameters. The error bars are 1 cr uncertainties, and do not 
include uncertainties from calibration 



faster in the outer parts of the disk (which are hardly detected in 
the weaker components) than in the inner parts. A close exam- 
ination of the shape of the surface, as a function of the ex- 
ponents p and indeed reveals for the strong hf components a 
broad, curved and relatively shallow global minimum, with two 
local minima that correspond to the two (flat and steep) solutions 
just described. 

We attempted fitting all components with a model including 
two power laws (crossing at some radius Rx), but did not find any 
sensibly better results: in most cases, the fit converged towards 
one of the two single power laws already identified. The rotation 
temperatures derived when varying R^ remained within the range 
of temperatures allowed by the single power-law fits as listed in 
TableSl 

Within the hypothesis where molecules are only in a sur- 
face layer, the depleti on parameter ckj can be constrained from 
the observations (see ISemenov et al.] |2Q08V the distribution of 
molecules being reasonably well described by two thin disks 
with slightly diflTerent inclinations, so that their projected kine- 
matic patterns differ. We find ad - ^ (c^d = ± 0.3 for the more 
sensitive CN measurements - see Table|4l), indicating that the de- 
tected molecules are most likely concentrated close to the disk 
plane. 

For the second approach, we must provide the underlying 
H2 density profile. We assumed two diflTerent disk models for 
the H2 density based on the analysis of the millimeter contin- 
uum emission performed by Guilloteau et al. ( 201 1). The first is 
a power law distribution; however, although the millimeter con- 
tinuum emission indicates outer radii of 200-300 AU, we ex- 
trapolate the H2 surface density up to 800 AU, which leads to a 
model with unlikely high densities at large radii. The second one 
uses the exponentially tapered solution of self- similar viscous 
evolution. Both solutions have similar surface densities near 100 
AU, but differ in the outer parts, i.e. beyond 200-250 AU. The 
H2 surface density is about 4 x 10^^ cm"^ near 100 AU. It de- 
creases as r"^-^ for the power law model, while the exponential 
taper leads to much steeper decline beyond 300 AU. For CN, the 
derived "depletion column density", is somewhat dependent 
on the disk model. For DlMTau, we find = 3^^^ x 10^^ cm"^ 
for the power law, and = 8 ± 1 x 10^^ cm"^ for the viscous 
model. In the former case, the error is highly asymmetric: the 3 
and 4 cr lower limits are : = 2 and 1x10^^ cm"^ respectively. 



For the Herbig star (IMWC 480), the sensitivity of our obser- 
vations is not good enough to derive both T and q. Nevertheless, 
assuming a slope q = 0, 3. value of T = 30 K is fitted for the CN 
emission. That is in accordance with the expected kinetic tem- 
perature of the molecular layer in the theoretical models. Using 
q = 0.4 instead results in a temperature r(300AU) = 14 K, 
but affects neither the derived surface density (by less than 5%), 
nor the exponent p (by less than 0.1). This is because in the 
10 - 30 K temperature range, the CN J=2-l l ine brightness is 
weakly dependent on the temperature (see Darto is et aP l2003l 
their Fig. 4, given for CO but also applicable to CN). For HCN, 
which has low signal-to-noise ratio, we assumed the same tem- 
perature (^ = 0, r = 30 K) and outer radius than derived from 
CN. Although arbitrary, this choice is justified by the results 
found for DlMTau and LkCal5, where CN and HCN display 
similar rotation temperatures and outer radii. 

The detection of HCN in MWC 480 is not easily illustrated, 
because the line intensity is spread over the three hyperfine com- 
ponents of the J=l-0 line. The best fit profile is shown in Fig.|4l 
Fig.[5] shows a signal-to -noise image obtain ed using the optimal 
filtering as described by lDutrey et al.l (l2007l) . i.e. weighting each 
channel by its modeled integrated intensity prior to averaging in 
velocity. There, the peak signal-to-noise is around 6, and clearly 
points toward the IMWC 480 star, thereby showing unequivocally 
the reality of the line emission from HCN. 

All temperatures may be increased by 10% to account for the 
calibration issues mentioned in Sec El An effect of similar mag- 
nitude, although not as linear as for the temperature, is expected 
on the molecular surface densities. However, the calibration un- 
certainties do not affect q, p, /?out, o;d or E^, which depend only 
the shape of the brightness distribution. 

4. Chemistry modeling 

Generic chemi cal models of cir c umstel lar disks have been 
published by "Aikawa & Herbsl^ ('1999'); "WiUacv & Langerl 
( 2000); Aikawa et al. (2002); van Zadelhoff et al. (2003); 
Aika wa & Nomural (120061) : IWiUacv et al.l (l2QQ6h : iFogel et all 
( 201 il). Each of these works intro duce a different degree of com- 
plexity. I Aikawa &Herbstl (11999!) consider gas phase chemistry, 
plus sticking and desorption, and assumed a vertically isother- 
mal disk model following the "Minimum Mass Solar Nebula" 
extended to 800 AU. The study from Willacv & Lanser ( 2000) 
uses the thermal structure derived from the . Chiang & GoldreictJ 
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some unknown properties of the observed disks (e.g. UV flux, 
mass, dust size), we performed new chemical network calcula- 
tions. The objective of these calculations is not to find an ap- 
propriate model for each of the observed disks, but to illus- 
trate common properties of disk chemical models, and pinpoint 
the dependency between the disk parameters and the predicted 
molecular column densities. We focus here on CN, which is pre- 
dicted to be formed mainly in the "warm" photo-dissociation 
layer, where chemical reactions are fast. This allows us to use a 
time-independent chemical model, rather than a more computer- 
intensive time dependent one, as would be needed for molecules 
formed deeper in the disk. 



Fig. 5. Signal to noise image for the HCN J=l-0 integrated emis- 
sion in MWC 480 (after removal of the continuum emission), 
with a velocity weighting profile equal to the best fit integrated 
model profile shown in Fig.|4l Contours are in steps of 2cr. 

Table 5. Elemental abundances with respect to total hydrogen 



4.1. FDR code 



Element 

C 

O 

N 

Mg 

S 

Si 

Fe 



abundance 
1.38 X 10-4 
3.02 X 10-4 
7.954 X 10-^ 
1.0 X 10-8 
2.0 X 10-6 
1.73 X 10-8 
1.7 X 10-9 



(1 19971) two-layer approxir nation, as well a s enha n ced photodes- 
orptio n yields following IWestley et al.l (Il995h . lAikawa et al.l 
(l2QQ2h used a D'Alessio disk model including vertical temper- 
ature gradients and self-consistently variable flaring, although 
dust and gas are assumed to be fully thermally coupled. All three 
models used a 1-l-lD approximation for radiative transfer, in 
which the stellar UV is attenuated along the line of sight to the 
star, while the ISRF impa cts isotropically on the disk surface. 
Ivan Zadelhoff' et al] (l2003l) improved the UV treatment, by using 
a 2-D radi ative transfe r code to solve for the UV field inside 
the disk, as'Foge l et al.l (|201 Shielding by H2 is treated in an 
approximate way, however, ^ikawa & Nomura ( 2006) further 
exp anded the na o dels b y considering the eff'ect of grain growth 
and ^Eogel et aP (l201ll) inves tigate the eff'ects of du st settling 
following th e prescription bvlD'Alessio et al.l (I2OO6I). W ith the 
exception of lFogel etaP (12011 *) and'WiUac v et al.l (12006). none 
of these models include any grain surface chemistry, which may 
be important in such environments. However, all of them predict 
that CN is constrained to a warm molecular layer, sometimes 
well above the disk plane, and that molecules are depleted in the 
disk plane. The column densities of CN and HCN do not vary 
much with radius. For a more direct comparison with with our 
results, we fitted the column densities predicted by the chemical 
models by power laws as function of radius, i.e. we derived 
and p. 

The model from ! Willacy et al.l (12006) is somewhat diff'er- 
ent, as it also includes turbulent diff'usion, and will be discussed 
later. A mor e elaborate set of g rain surface chemistry has been 
included by IWalsh et aP (12010 ^). who also find similar layered 
structure. 

Existing chemical models are not necessarily tailored to the 
sources we have observed. To better understand the impact of 



We use the PDR code from the Meudon group (iLe Bourlot et al.l 
ll993l:lLe Petit et al.ll2006l) n iodified to account for a non uniform 
grain size distribution (see IChapillon et 10 120081) . The chemi- 
cal model assumes a one-dimensional stationary plane-parallel 
slab of gas and dust illuminated by an ultraviolet (UV) radiation 
field. The radiative transfer in the UV takes into account self 
and mutual shielding of H, H2 and CO lines and absorption and 
diff'usion of the continuum radiation by dust grains in a 1-D ge- 
ometry. Molecular abundances and, optionally, thermal balance 
of the dust grains and the gas are calculated iteratively at each 
point in the cloud. The chemical network is similar to that of 
iGoicoechea et al.l (l2006l) . The chemical model does not consider 
the freeze-out of atoms and molecules onto grains and surface 
chemistry reactions, to the notable exception of those leading to 
the formation of H2. 

The assumed grain size distribution is a standard power-law 
n{a) oc a"^, where a+ and a- are the maximum and minimum 
cutoff' radii, respectively. The introduction of a range of grain 
sizes aff'ects the UV extinction curve, the rate of H2 formation 
and the thermal balance of the gas. The resulting extinction curve 
is calculated using the Mie theory for homogeneous isotropic 
spherical particles. The composition of the dust is 50% silicate 
and 50% graphite. We assume y = -3.5 (the value observed 
in nearby interstellar clouds) and a- = 3nm in all calculations. 
This minimum grain radius a- = 3nm is chosen small enough to 
properly account for the photoelectric heating process, the UV 
extinction, and the formation of H2. 

Keeping the dust mass constant, we have simulated grain 
growth by varying the maximum radius a+(03/dm - 0.1mm) 
while keeping the exponent y constant for the sake of simplicity. 
The amount of small grains therefore is reduced to the benefit of 
large grains, modifying accordingly the extinction curve. With 
y = -3.5 kept const ant, the UV opa city scales as 1/ ^/a^ for 
a+ > lO^m (^see IChapillon et al.ll2008l) . 

The chemistry code is one-dimensional and calculates the 
temperature, radiation field and molecular abundances in the di- 
rection perpendicular to the surface, that we assume here perpen- 
dicular to the disk plane. The justification for such an approach is 
that small dust grains in the upper disk atmosphere scatter (half 
of) the incident UV ffux towards the disk. We mimic the second 
dimension by resuming the calculation at diff'erent disk radii (12 
to 15 diff'erent radii from 50 AU to the external radius Rput), the 
output being a 1-hlD model (^ see IChapillon et al]|2008l for more 
details). 
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4.2. Set of PDR models 

We explore a range of parameters to understand the effects of 
UV field intensity, grain size and cosmic ray ionization rate. 
Three basic disk models (with diff'erent density and temperature 
structure) are used to represent the 3 sources DMTau, LkCa 15 
and MCW480. The in put disk structure is similar to those of 
ID' Alessio et alJ (Il999h . The density and (dust) temperature dis- 
tributions were computed self-consistently from the disk masses 
and stellar properties of each source by D.Wiebe & D.Semenov 
(private communication) assuming uniform grain size of 0.1 yum; 
the disk models display a vertical te mperature gradient, the tem - 
perature increasing with height (e.g. lChiang & Goldreichll 1997 1). 
Except when otherwise noted, we assumed the gas temperature 
to be equal to the dust temperature. The elemental abundances 
used are given in Table [5] 

For each disk, we vary several parameters. 

Grain distribution: We have simulated grain growth by 
varying the maximum radius a+\ in the "big grains" case 
a+=0.1 mm, in the "small grains" case a+=03 yum. 

The incident UV field Illuminating the disk surface is the 
sum of the interstellar UV field (ISRF, due mainly to the sur- 
rounding massive stars) and the stellar UV field coming from 
the central star. T-Tauri spectra present strong UV excess due 
to shock accretion of material onto the star. This UV field is 
assumed to follow the shape of an ISRF Draine field with an in- 
tensity scaling factor;^ at r = 100 AU. The UV field decreases 
with in creasing radius as 1/r^. Acco rding toiHerbig & GoodrichI 
(Il986h X = 10"^ is appropriate, but iBergin et al.l (l2003h find a 
much lower value, y = 500. Using the measurements from 
iBergin et al.l (l2003b . iHenning et al.l (l2010h derived UV fluxes 
for DMTau (x = 300), LkCal5 (x = 1850) and MWC480 
(X = 5000) (see their Table 1). Accordingly we explore the range 
X = 10^-10^ and run also some additional models with;^ =10^ 
to see the impact of the stellar UV field intensity on the CN/HCN 
chemistry. 

In addition, a couple of models were run for DM Tau with 
the photodissociation rates calculated with out UV excess, but 
for a b lack body at 4000 K using rates from Ivan Dishoeck et al.l 
(12006). Results are presented only in Fig.[8l 

The Cosmic Ray ionisation rate (CR): We performed cal- 
culations with the standard value of 10"^^s"^ and with an in- 
creased rate of 5 10"^^s."^. 

Physical structure: In addition to the values intrinsic to the 
basic disk models, we vary the density (by a factor of 3) and the 
temperature (by ±30%) around those values. In the later case, 
the scale height is not recomputed, but kept equal to that of the 
basic model. 

The thermal balance was computed only at R=300 AU in a 
few models (those presented on Fig.[TO|, the incident UV field 
being the main ener gy source. Th e dust temperature was de- 
rived following Burt on et al.l (11990 ) and used to compute the 
energy exchange betwe en the gas and the grains, as described 
in Burke & Hollenbact] (Il983l) . The gas temperature is then cal- 
culated, taking into account heating and cooling processes (see 
iLe Petit et al.ll2006l for more details). Again, the scale height is 
that of the bas ic model and is not modified consistently. 

Ly a line: IBergin et^ (l2003h pointed out the possible im- 
portance of the Ly a line on the chemistry. Essentially, HCN can 
be photo-dissociated by Ly a radiation while CN cannot. As the 
Ly a line can be very strong in the spectra of T-Tauri (carrying 
up to a few 10% of the UV energy), this aff'ects the photodissoci- 
ation rate, and thus the CN/HCN chemistry. To test the influence 
of Ly a radiation, we implemented its impact on the chemistry 



following Ivan Dishoeck et al.l (l2006l) . and ran a few models in 
which the Ly a line represents l/3rd of the total UV luminosity. 

Unless otherwise specified, the model parameter values are 
as follows: a+=03jim (i.e. small grains), a ISRF UV field with 
;^ = 10^ and no hya included, CR=10"^^s"\ nominal density 
and temperature and no thermal balance calculated. 

4.3. PDR modeling results 

The column densities predicted by the chemical model are fitted 
by power laws as functions of radius, and are thus characterized 
by the values of the S300 and p, like for the observations. Power 
laws off'er fairly good approximations to the resulting distribu- 
tion of column densities. 

The calculated S300 and as a function of the various input 
parameters are illustrated for DM Tau in Fig. [6] for CN and Fig.|7] 
for HCN. The impact of the model parameters is shown in a 
more compact form for all sources in Fig. IB. II A comparison 
with other chemical models, fitted in the same way, is displayed 
in Fig. [8] together with the observed values. 

The main results are: 

- The calculated column densities (at 300 AU) E300 match rea- 
sonably well the observed values (within a factor of a few for the 
best set of parameters), but the radial dependence (see parameter 
p) is weaker in the models (p(mod) < p(obs)). 

- The most influential parameters are the UV flux and the grain 
sizes. These two parameters control the photodissociation pro- 
cess. 

- As a consequence a high CN/HCN ratio is observed for large 
grains. 

- A modification of the CR influences the HCN/CN chemistry 
only in the case of small grains, i.e. when UV penetration is ef- 
ficiently blocked. 

- Despite the lack of molecule freeze out and grain chemistry, our 
PDR code results for CN and HCN are in broad agreement (flat 
radial distribution and similar value of column densities within a 
factor of a few) with publis hed results of models that take these 
processes into account (see Aikawa & Herbst (1999, Fig. 7-9), 
FAikawaetal. (200 2, Fig^ 3), van Zadelhofl' et al. (2003, Fig. 5 
and 8). I Aikawa & Nomural (12006 . Fig. 9)). This suggests that the 
freeze-out of molecules onto grains, as implemented in these 
works, does not afl'ect significantly HCN and especially CN, 
which remain domina ted by ph oto-chemistry in these models. 
Note, however, that Walsh et al.l ( 2010) find out the grain surface 
chemistry decreases the HCN column density, but only mildly 
aff'ect CN. 

As in the other chemical studies, the CN molecule is located 
within a narrow layer close to the surface (Ay -0.1 correspond- 
ing to z~ 200 AU (n~ 5 x 10^ cm"^) in the small grain case 
and z~ 50 AU (n~ 7 x 10^ cm"^) in the big grain case) whereas 
HCN is buried more deeply inside the disk. The vertical distri- 
bution of molecules at r = 300 AU are presented in Fig. lB.2l for 
small and big grains and two diff'erent values of the incident UV 
field. As we do not include sticking onto grains, HCN remains 
fairly abundant in the mid-plane contrary to models including 
sticking onto grains where HCN is not present under z/r -0.1 
(e.g., Aikawa & Nomura 2006). 

Relative abundances can be a misleading indicator of the lo- 
cation of the bulk of molecules, because of the strong vertical 
density gradient. To better indicate where molecules appear, we 
trace in Fig. [9] the cumulative column density from the disk mid- 
plane to the surface as a function of the gas density. The latter 
can be converted into height above the disk mid-plane using the 
vertical distribution. The temperature (which, in this model, ap- 
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Fig. 6. Variations of the radial distribution of CN ( S300 and p parameters) according to the maximum grain size a+, the UV field 
intensity ix factor), the rate of cosmic ray ionization and modification of the physical structure. The "standard" values are;^ = 10^, 
CR=10"^^s"^ and nominal temperature and density structure. Squares: big grains ((2+=0.1 mm), triangles: small grains ((2+=0.3 yum). 



plies equally to gas and dust) is also displayed. Fig. [9] shows the 
temperatures and densities at which the bulk of the column den- 
sity is built up, allowing to predict what is the average line ex- 
citation temperature. For essentially all models the CN column 
density is mostly build up in a region where the kinetic temper- 
ature is > 30 K and at a density n > lO^cm"-^, the only exception 
being the small grains, low UV (x = 10^) case, where only half 
of the CN column density is in these conditions. 

Similar results are obtained in models for which the gas and 
dust thermal balance is computed: the vertical distributions of 
CN and HCN abundances are not significantly aff'ected, although 
the gas temperature is in general somewhat higher in this case 
(Fig.[TO|). The chemistry is not significantly aff'ected because 
most reactions depend weakly on the temperature in this tem- 
perature range. 

5. Discussion 

The results from our power law model fitting of the observations 
of the T-Tauri stars DM Tau and LkCal5 appear to contradict 
the chemical model predictions on several aspects. 

1 . While ^^CO indicates temperatures in the rang e 20 to 30 K at 
100 AU for the upper layers (.Pietu et al.ll2007.) . the derived 



CN and HCN rotation temperatures, which apply mainly to 
radii 100-300 AU, where our observations are the most sen- 
sitive, are very low (8 to 12 K, including calibration errors). 

2. In addition, the constraint on ad - (equivalent to large Sj) 
suggests that CN molecules appear closer to the disk than ex- 
pected. This is coherent with a low temperature but in appar- 
ent contradiction with the location of the predicted molecular 
layers. 

3. Most models predict slopes p which are much lower 
than observed; the ma in exception being the results of 
lAikawa & HerbstI (Il999l) for CN. 



For the Herbig star (MWC 480), the temperature may be ade- 
quate, but the slope problem remains. This leads to several open 
questions : 



1 . Do we understand the basic physical structure (density, tem- 
perature) of protoplanetary disks correctly ? 

2. How a misleading structure can aff'ect the line formation ? 

3. Do we understand the chemistry, in particular that of CN 
well enough ? 
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Fig. 7. Variations of the radial distribution of HCN (S300 and p parameters) according to the maximum grain size a+, the UV field 
intensity (x factor), the rate of cosmic ray ionization and modification of the physical structure. The "standard" values are;^ = 10^, 
CR=10"^^s"^ and nominal temperature and density structure. Squares: big grains (a+=0.1 mm), triangles: small grains ((2+=0.3 yum). 



5A. Excitation conditions 

CN is not the only photo-dissociation sensitive molecule dis- 
playing surprising low excitation tempe rature s . From a spatially 
resolved, two-transition study, Henning et alj ( |2Q1Q| ) found that 
C2H display the very same behavior, which cannot be repre- 
sented with current chemical models. While for CN the debate 
remains possible, the chemistry of C2H is relatively well known, 
so the coincidence points towards a general problem to retain 
molecules at low temperatures in proto-planetary disks. 

In all chemical models, HCN is formed closer to the mid- 
plane than CN, at densities largely sufficient to thermalize 
the HCN J=l-Q transition (5 x 10^ cm"^ G uilloteau & Baudrv l 
(Il98lh : iMonteiro & Stutzkil (Il986h ) The HCN excitation tem- 
perature derived from our observations is then expected to re- 
flect the kinetic temperature of the disk mid-plane in the region 
where the J= 1-0 line is optically thick, that is up to about 150 AU 
in DM Tau and LkCa 15. The values derived from the observa- 
tions are ^ 6 K for DM Tau and ^ 12 K for LkCa 15, the excita- 
tion temperature for MWC 480 remaining unconstrained. CN, on 
the other hand, is expected to be mainly formed (at least largely 
formed) in the upper layer, close to the disk surface. Although 
the overall finding that r(CN) > r(HCN) in both sources is in 



qualitative agreement with this expectation, the observed values 
are surprisingly low. 

From Fig. [21 for all models with small grains, 80 % of the CN 
column density builds up in a region where the kinetic tempera- 
ture is > 30 K. Furthermore, 50 % of this column density origi- 
nates from regions where the density is at least 5 - 8 x lO^cm"^, 
and CN is not present at all for densities exceeding about 
3 X lO^cm"^, i.e. below 2 scale heights. This is in contradiction 
with the observations which indicate that the bulk of CN is near 
the disk plane (ad < 1 at the 3 cr level). 

While all published disk chemical models clearly have high 
CN abundances only at large heights, deriving the densities and 
temperatures at which the bulk of the CN column density is built 
is often not possible, as most studies focussed more on the gen- 
eral behavior of the chemistry than on detailed l i ne em ission pre- 
diction. Nevertheless, we note that lWalsh et al.l (l2010h only have 
substantial amount of CN for T > 20 K and n > 10^ cm"^, at 
rather large z/r, as in our model. 

The comparison with other transitions, CN(l-O) and CN(3-2) 
show that for MWC 480, our model predict a stronger CN(3-2) 
line (9.0 Jy.km/s) than observed. This can be ascribed to partial 
sub-thermal excitation for this transition, as its critical density 
is 3 - 6 X 10^ cm"^. For DM Tau, our predictions for the CN 
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Fig. 9. Cumulative CN (red) and HCN (blue) column densities 
and temperature (black) from the mid-plane to the atmosphere 
as a function of the density at a radius of 304 AU in the DM Tau 
disk. Models with temperature fixed from the structure file. This 
samples only half of the disk, so the column densities have to be 
multiplied by 2 to retrieve the S values. 
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Fig. 10. Same as Fig. [9] for models with thermal balance com- 
puted. 



temperatures below 15 K by significant sub-thermal excitation 
appears difficult. 



J=l-0 agree with the observations of lDutrev et aP (Il997l) . but 
there is no published data on the CN(3-2) line. For the disk of 
LkCal5, our best fit model underestimates by 30 % the line flux 
of the CN(3-2) line. Raising the rotation temperature to about 
25 K brings the predicted line flux in better agreement for this 
transition. This suggests that a fraction of the CN molecules lie 
in a warmer (and dense enough) temperature region. However, 
although the exact collision rates for CN with H2 are unknown, 
the densities are high enough to be close to thermalization for 
the observed lines, and it is likely that T^ot > Tk/2, so explaining 



5.2. Surface densities of CN and HCN 

A general prediction of all chemical models is that the surface 
densities of CN and HCN are nearly constant as a function of 
radius p ^ 0. The observed slopes are larger, ^ 1, and the sep- 
arate hyperfine CN analysis even suggests a steepening of the 
surface density distribution at larger radii. Furthermore, chang- 
ing the model disk mass to 1/3^^ or 3 times the nominal value 
only changes the surface density of CN and HCN by about 10- 
20 % (see Fig.[6l[7landlBJl). 
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Fig. 11. CN/HCN column density ratios as a function of the UV 
field at R = 97AU (triangles), 304 AU (squares), and 602 AU 
(circles) for the models with small grains (thick black lines) and 
big grains (thin red lines) 



iBergin et"aD (l2003l) showed that Ly a radiation significantly 
aff'ects the relative photodissociation rates of HCN and CN. In 
our model, the impact of the Ly a line on CN is quite small. 
This is because the direct production of CN from the photodis- 
sociation of HCN (HCN + /zv ^ CN + H) is not in general 
the main production path of CN. Figure [121 presents the main 
chemical reaction leading to the formation of CN at a radius of 
304 AU. Photodissociation of HCN is the most efficient CN pro- 
duction mechanism only in a small region with high HCN abun- 
dance. Upwards in the disk, the main CN formation reactions are 
N CH ^ CN H and N -h C2 ^ CN -h C, so that the over- 
all abundance of CN is controlled by the abundance of atomic 
Nitrogen. However, HCN decreases with Lya intensity, so that 
the CN/HCN rat io depends on the Lya flux (as also found by 
iFogel et aDl2011i) . We find however that the effect is very small 
when dust is composed of small dust grains. With large grains, 
the HCN column density is divided by a factor ^ 2.5 for a Ly a 
line contributing to l/3rd of the UV luminosity. This is consistent 
with the lFogel et al.l (|201 1) results where the effect is of Ly a is 
visible only for settled grain models (i.e. models with a reduced 
UV extinction in the upper layers, as for large grains). 

Finally, PDR models also predict that the CN/HCN ratio in- 
creases with grain size because large grains allow UV flux to 
penetrate deeper in the disk. This effect is illustrated in Fig.fTT] 
The high CN/HCN ratio (4 to 10 at 300 AU) in the three sources 
is in agreement with the presence of large grains in the three 
disks. On the other hand, Aika wa & Nomural (120061) results sug- 
gest a rather weak dependency of the CN/HCN ratio on the grain 
sizes. However, this study explored much bigger grain sizes than 
we do here, but the differences between their 10 fim and 1 mm 
case is similar to our findings. Walsh et al. ( 2010) also found that 
the CN/HCN ratio increases when surface chemistry is included, 
as Nitrogen can be incorporated into larger molecules in the disk 
mid-plane. 

5.3. Toward a cold molecular Layer ? 

We investigate here mechanisms which would partially provide 
cold molecular material close to the mid-plane. 
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Fig. 12. Efficiency of CN formation reactions at R=304AU 
for DM Tau with smafl grains, ;^ = 10^ CR=10-^^s-^ Only 
the main reactions are plotted (i.e. more than 10%). 1 (red): 
H -h CN+ ^ H+ -h CN, 2 (green): N -h C2 ^ CN -h C, 3 (blue): 
N -h CH ^ CN -h H, 4 (cyan): HCN + photon ^ CN -h H, 5 
(yellow): HNC -h photon ^ CN -h H, 6 (magenta): OCN -h pho- 
ton O -hCN. The cumulative CN column density is plotted in 
black. 



5.3.1 . Large grains 

Grain growth can bring molecules closer to the disk mid-plane, 
by allowing UV radiation to penetrate deeper in the disk struc- 
ture. Fig. [9] also indicates that the only case in which CN is pro- 
duced in sufficient quantities at low temperatures is the high UV 
ix = 10"^), big grains a+ =0.1 mm case, where 50% of the CN 
column density builds up at temperatures below 10 K, and the 
rest between 10 and 20 K. However, neglecting freeze-out and 
surface chemistry in this case is unlikely to be valid. Moreover, 
in this case, the disk temperature and density structure was im- 
posed a priori, and is not consistent with the impinging UV field 
and its attenuation. If, using the same density structure, we com- 
pute the thermal balance assuming the incident UV field is the 
only heating source, the picture changes completely (Fig.fTOl). 
The whole gaseous disk becomes warm, the minimum tempera- 
ture being 60 K. Furthermore, 90 % of the CN column density is 
build up at densities above 10^ cm"^, thus sub-thermal excitation 
becomes also unlikely. 

Enhanced photo-desorption would bring larger amounts 
of (cold) CN, H CN and C2H near the disk mid-plane. 
IWillacv & LangeJ (l2000h checked this possibility using moder- 
ate grain growth but it does not provide enough cold gas. Doing 
the same calculation within larger grains appears however as an 
interesting clue to partially reconcile observations and models. 



5.3.2. Turbulent mixing 

Turbulent transport could bring CN and C2H formed in the 
warm upper layers down to the disk mid-plane, where they 
would cool very quickly. The main difficulty is to avoid stick- 
ing onto cold dust grains. Hersant et al. ( 2009) explored this for 
CO molecules, and found that turbulent diffusion does not affect 
much the total amount o f CO, but cou l d brin g a small fraction 
of it to low temperatures. I Willacy et al.l (l2006h managed to get a 
substantial amount of CN towards the disk plane, but only using 
a large diff'usion coefficient, 10"^^ cm^.s"\ which is 10 times 
larger than expected from our current estimates of the viscos- 
ity (a viscosity parameter below 0.01) in disks. We note how- 
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ever that the model predictions are (at least for CN) not too 
far from the observed column densities. The low mixing effi- 
ciency, which is a result of dilution (the gas density in the upper 
laye r being small in com parison with the disk mid-plane, see 
also lSemenov et al" may thus not be a major problem, as 

we only want to bring the existing CN molecules to low temper- 
ature regions, but not to increase significantly the total amount 
of CN. 

5.3.3. Lowering the gas-to-dust ratio 

In our resolved observations, is actually geometrically con- 
strained by measuring the height above the disk mid-plane where 
the molecules are located; it is expressed as a column density 
by simply scaling to the total disk surface density (see Eql6]). 
When assuming A:y(230GHz) = 2 cm^.g"^ (of dust) and a canon- 
ical gas-to-dust ratio of 100, the measured "depletion" column 
density, Sj, appears -4-10 times larger than predicted from 
chemical models (see in particular Aikawa & Nomura 2006). 
[Aikawa & Nomural (|2006) also show that does not strongly 
depend on the grain size: so what controls the location of the 
CN layer (and those of other molecules) is essentially the col- 
umn density of hydrogen towards the disk surface. Lowering the 
gas-to-dust ratio to 20 (or reducing the hydrogen surface density 
by a factor 5), we bring to first order the observed in rough 
agreement with the model prediction. In such a case, CN would 
be present in regions where the density is low enough (a few 10^ 
cm"^) to lead to substantial sub-thermal excitation. 

We checked this possibility by using our PDR and radiative 
transfer codes in a disk 6 times less massive than our nominal 
assumption. The excitation is computed using the modified LVG 
approximation as described in Pavlyuchenkov et al. ( 2007) who 
also show that this is a reasonably accurate solution for the radia- 
tive transfer in disks. The collision r ates for CN were assume d 
identical to those of HCN with He (iGreen & Thaddeusl[T974h . 
For DM Tau, using the tapered edge surface density profile, we 
indeed find out = 1. 7 + 0. 1x10^^ cm"^, and a kinetic temper- 
ature of 15 - 20 K. Given the uncertainties on the collision rates, 
such a solution appears compatible with the expected tempera- 
ture range in the molecular rich region. Such a low gas-to-dust 
ratio would also provide an appropriate explanation fo r the l ow 
temperatures also derived from C2H by Henn ing et al.1 (l2010l) . 

Unfortunately, with our chemical model, such low disk 
masses result in predicted HCN and CN column densities which 
are substantially lower than observed, so that this mechanism is 
not fully satisfactory to explain the observations. Imaging higher 
lying transitions, which have higher critical densities than the 
J=2-l line observed here, could give more accurate constraints. 

Note that, rather than changing the gas-to-dust ratio, rais- 
ing the mm wave dust absorption coefficient by a similar factor 
is also a strictly equivalent alternative, and leads to the same 
low disk masses (0.003 M© for DM Tau). However, our assumed 
value of /c(1.3mm) = 2cm^.g"^ (of dust) already appears large 
(see e.g. lDraindl2006, and references therein). 

5.3.4. Low temperature chemistry 

Finallv. iHifv-Blant et al] (l2008h have observed very cold ^^CN 
in the nuclei of a few dense cores of known temperature and 
density (T ~ 6 - lOK, /t ~ 10^ - lO^cm"^), which present 
temperature and densities relatively similar to those encoun- 
tered in the mid-planes of TTauri disks. They find that CN is 
not strongly depleted in the nuclei of the cores, contrary to 



CO, and can be a fai r column density tracer in these extreme 
conditions. However, Hily-Blant et al.' ("2010) were unable to 
account for the near constancy of the CN/HCN ratio in their 
chemical models. These observations lead also open the pos- 
sibil ity that cold chemistry is not fully understood at least for 
CN. iHilv-Blant et al.1 (l2010l) conclude that our knowledge of the 
low temperature chemistry of CN is limited by the extrapola- 
tion of the rate coefficients measured at high temperatures to 
the 10 K regime for many of the important reactions, includ- 
ing those mentioned in S ec l5.2l (iPineau des Forets et al.l I199QI: 
iBoger & Sternbergll2005l) . However, at such low temperatures 
and high densities, chemist ry also depe nds on reactions on the 
grain surfaces, as shown by I Walsh et al.l (|2010) in particular for 
HCN. Hence, another viable alternative would be inappropriate 
knowledge of the grain surface chemistry. 

6. Summary 

We present sensitive, high spatial and spectral resolution obser- 
vations of HCN J=l-0 and CN J=2-l in 2 T Tauri disks (DM Tau 
and LkCa 15) and one Herbig Ae (MWC 480). Column densities 
and excitation conditions are recovered through a minimization 
scheme which utilizes all the information provided by the hyper- 
fine structure of the observed transitions. Although we find out 
r(CN) > r(HCN) as expected from a layered disk structure, the 
analysis indicates that the location of CN and its apparent ex- 
citation temperature appears incompatible with the expectations 
from current chemical models. Lowering the gas-to-dust ratio 
might partially solve the conflict on temperature and localiza- 
tion of the molecules, but may not produce large enough quanti- 
ties of CN and HCN. The observations thus suggest that a sub- 
stantial fraction of the molecules are created in the disk plane at 
low temperatures, perhaps by surface chemistry. Sensitive, spa- 
tially resolved images of higher excitation lines of CN and HCN 
would be required to provide a complete diagnostic to determine 
the excitation conditions in disks allowing us to decide between 
the possible options. 
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Fig. A.l. Channel maps of CN hyperfine components towards DM Tau. From top to bottom, the velocity scales are referred to 
226.874745 GHz (3 blended components), 226.679382 GHz, 226.663703 GHz and 226.659575 GHz. Note the better spectral 
resolution in the top panel, the spatial resolution is 1.7 x 1.2'' at PA 120°, contour spacing is 53 mJy/beam, or 0.61 K, corresponding 
to 2. Oct in the upper panel, and 23cr in the lower ones. 
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Mg. A.2. As Fig jA.U but towards LkCal5. All panels use the same spectral resolution. Note the repeated channels in the two bottom 
panels, as the velocity spacing between these components is 5 km/s, i.e. 12 channels. The spatial resolution is 1.7 x 1.0'' at PA 44°, 
contour spacing is 60 mJy/beam, or 0.81 K and 2.lcr. 
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Fig. A.3. As Fig jA.2l but towards MWC 480. The spatial resolution is 2.2 x 1.35'' at PA 60°, contour spacing is 55 mJy/beam, or 
0.45 K and 2cr. The large rotation velocity result in an overlap of emission between the two hyperfine components in the bottom 
panels. 
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Fig. A.4. Channel maps of HCN hyperfine components towards DMTau. The spatial resolution is 3.7 x 2.4'' at PA 39°, contour 
spacing is 12 mJy/beam, or 0.21 K, 2. Oct. 
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Fig. A.5. Channel maps of HCN hyperfine components towards LkCal5. The spatial resolution is 4.1 x 2.9'' at PA 51°, contour 
spacing is 16 mJy/beam, or 0.21 K and 2. Oct. 
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Fig. A.6. Channel maps of HCN hyperfine components towards MWC480. The spatial resolution is 4.5 x 3.6'' at PA 66°, contour 
spacing is 16 mJy/beam, or 0.15 K and l.Scr. 
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Fig.B.l. Value of X300 and p from our chemistry modeling for the three sources, with our "standard" chemistry, o: density and 
temperature from theoretical models, p: temperature increased by 30%, m: temperature decreased by 30%, f : density multiplied by 
a factor 3 and d: density divided by a factor 3. In many models the blue and green marker overlap, as well as the red and black ones 
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Fig. B.2. H, H2, C+, C, CO, CN, HCN abundance and CN/HCN ratio as a function of opacity from atmosphere (Av=0) to mid-plane 
at R-304 AU for DMTau with CR^IQ-^^ s'^. 
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